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Introduction 


In the past, the design of rocket engines has primarily relied on the cold flow/hot fire test, 
and the empirical correlations developed based on the database from previous designs. 

However, it is very costly to fabricate and test various hardware designs during the design cycle, 
whereas the empirical model becomes unreliable in designing the advanced rocket engine where 
its operating conditions exceed the range of the database. The main goal of the 2 nd Generation 
Reusable Launching Vehicle (GEN-II RLV) is to reduce the cost per payload and to extend the 
life of the hardware, which poses a great challenge to the rocket engine design. Hence, 
understanding the flow characteristics in each engine components is thus critical to the engine 
design. In the last few decades, the methodology of computational fluid dynamics (CFD) has 
been advanced to be a mature tool of analyzing various engine components. Therefore, it is 
important for the CFD design tool to be able to properly simulate the hot flow environment near 
the liquid injector, and thus to accurately predict the heat load to the injector faceplate. 

However, to date it is still not feasible to conduct CFD simulations of the detailed flowfield with 
very complicated geometries such as fluid flow and heat transfer in an injector assembly and 
through a porous plate, which requires gigantic computer memories and power to resolve the 
detailed geometry. The rigimesh (a sintered metal material), utilized to reduce the heat load to 
the faceplate, is one of the design concepts for the injector faceplate of the GEN-II RLV. In 
addition, the injector assembly is designed to distribute propellants into the combustion chamber 
of the liquid rocket engine. A porosity mode thus becomes a necessity for the CFD code in order 
to efficiently simulate the flow and heat transfer in these porous media, and maintain good 
accuracy in describing the flow fields. Currently, the FDNS (Finite Difference Navier-S takes) 
code is one of the CFD codes which are most widely used by research engineers at NASA 
Marshall Space Flight Center (MSFC) to simulate various flow problems related to rocket 
engines. The objective of this research work during the 10-week summer faculty fellowship 
program was to 1) debug the framework of the porosity model in the current FDNS code, and 2) 
validate the porosity model by simulating flows through various porous media such as tube 
banks and porous plate. 

Numerical Methodology 

The framework of the FDNS code is an elliptic, finite difference Navier-Stokes flow 
solver " , which employs a predictor plus multi-corrector pressure-based solution algorithm. 
Higher order upwind, total variation diminishing (TVD), or central difference schemes plus 
adaptive second-order and fourth-order dissipation terms are used to approximate the convection 
terms of the transport equations. Various matrix solvers, such as vectorized point implicit, 
conjugate gradient, and generalized minimal residual 4 (GMRES), are provided in the code such 
that users can select one for a given transport equation. Finite -rate and equilibrium chemistries 
were incorporated in the code to account for the effect of combustion. Since the FDNS flow 
solver is a structured code, multi-block, multi-zone options are included in the code so that 
problems with complex geometries can be analyzed efficiently. The governing equations in 
curvilinear coordinates solved in the code can be expressed in general form as 

1 dpg _ dF t | s 

J dt q 
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where q stands for dependent variables (1, u, v, w, h, k, s, and a n which denote unity, velocities 
in the Cartesian x-, y-, z-direction, enthalpy, turbulence kinetic energy and dissipation rate, and 
species mass fractions respectively). £/, U and p are the curvilinear coordinates, time and fluid 
density. The numerical flux in the i-direction, F { , consists of convective flux, Fc„ and a viscous 
flux, Fvu i.e. 

F i = Fc j +Fv i , F Cj = pU jq , Fv i = -p e G ij ^~ 

where U u Gy, and / represent the transfonned velocities, diffusion metrics, and Jacobian of the 
coordinate transformation, respectively. The source terms, S q , of the continuity, x-, y-, z- 
momentums, energy, turbulent kinetic energy and dissipation rate, and species equations 


Porosity Model 


The porosity model is employed to account for the effects of area/volume blockage and 
drag force as well as heat source/sink due to the presence of the porous media without resolving 
the detailed geometry of tiny pores and solid objects. The implementation of the porosity into 
the CFD code is illustrated as follows. First, the volume porosity (y v defined as the ratio of the 
volume occupied by the fluid to the total volume), and the surface porosity (y defined as the ratio 
of the surface area in the i-direction available for the flow passage to the corresponding total 
surface area in the same direction) are introduced into the governing equation and can be 
expressed as 


Yv dpq d Yi F i 


+ y v S .. 

J dt q 

Second, a distributed drag force, D, and a heat flux source/sink tenn, //, are added to the 
right hand side of the momentum equation and the energy equation, respectively, to account for 
the effect of resistance and heat transfer due to the presence of the porous media. These drag 
force and heat flux terms were modeled based on geometric parameters and the averaged 
velocity through the porous media. The drag force D is defined as 


d = \ p v; 


c d a 


where Cd is the local drag coefficient is a function of types of porous media, local flow speed 
and geometric parameter. For the flow through an inline or staggered tube bank, empirical 
correlations 5 in Tables 1 and 2, are incorporated into the code to compute the drag coefficient. 


Table 1 : Empirical correlations of the drag coefficient for the flow through in-line tube banks 

Range 

C D 

10 2 < Re</< 4 x 10 3 

0.535 exp(5.378 Re/' 345 ) 

4x 10 3 <Re rf <6x 10 4 

0.735 -0.411 x 10' 6 Red 

6x 10 4 <Re rf < 1 x 10 6 

0.621 + 0.169 x 10' 5 Re^- 11.343 x 10' 12 Re/ + 16.656 x 10' 18 Re/ 
-7.377 x 10' 24 Re/ 

1 x 10 6 <Re rf 

0.2735 
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Table 2: Empirical correlations of the drag coefficient for flow through staggered tube banks 

Range 

C D 

10 2 <Re rf <4 x 10 3 

0.417 exp(4.932 Re/ 296 ) 

4 x 10 3 < Re^< 6 x 10 4 

0.647- 0.5 x 10' 6 Re</ 

6 x 1 0 4 < Rej < 1 x 10 6 

0.618 + 0.491 x 10' 6 Rerf- 6.303 x 10' 12 Re/ + 10.694 x 10' 18 Re/ 
-5.2 x 10' 24 Re/ 

1 x 10 6 < Re^ 

0.2735 


Both the drag coefficient model for a perforated plate and the heat flux model are incorporated 
into the code. These models are detailed in the final report to NASA MSFC TD-64. Due to the 
space limitation of this report, these models will not be described here. Readers can refer to the 
final report to TD-64 for the detailed. 

Numerical Results 

Tube Banks 


Three flow conditions (Re</ = 1000, 5000, and 10000) were simulated for the flow 
through both in-line and staggered tube banks with three different tube spacing ( a = Si Id = 2, 1.5, 
and 1.25), where the tube has a diameter of 1 cm. Due to the space limitation of this report, only 
the result of the staggered tube bank flow is included here. The complete comparison and 
discussion are detailed in the final report to TD-64. The computational domain and definitions 
of key parameters are sketched as shown in Figure 1. The comparisons of the Euler number 

(Eu = 2 A p / pu 2 ) between the numerical results and the empirical data are plotted as shown in 
Figure 2. Though the pressure loss is shown to be dependent on the flow blockage; it does not 
demonstrate the similar trend for various Reynolds number. Later, it was found out the empirical 
correlation of drag coefficients has less dependence on the Reynolds number than the pressure 
drop data. Hence the empirical correlation of drag coefficients needs to be fine-tuned in the 
future to obtain the Reynolds number dependence. 



parameters of the staggered tube bank 


Pressure Loss for the Flow Through Staggered Tube Bank 



Figure 2: Pressure loss coefficient for the 
flow through staggered tube banks 
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Porous Plates 


Numerical simulations were conducted for the flow through a l”-thickness porous plate 
(pore diameter = 10 pin, volume porosity = 0.1) with various specific mass fluxes. The 
numerical result indicates the pressure loss predicted by the porosity model demonstrates a trend 
similar to the test data for various specific mass fluxes. Due to the proprietary of the test case as 
well as the space limitation of this report, the data comparison will not be presented here. If 
interested, readers can refer to the final report to NASA MSFC TD-64. 

Conclusion 


The framework of the porosity model in the FDNS-500 code has been verified. The 
numerical analyses of the tube bank flow demonstrate the qualitative trend of the pressure loss 
across the tube bank. However, more validations are needed to fine-tune the model to predict 
both qualitative and quantitative trends for various area porosities and Reynolds numbers. More 
validations are also needed for the porosity model of the perforated plate to examine the effect of 
different pore sizes and volume porosities. Due to the limitation of the time, only the empirical 
correlation for the drag force of the porosity model was verified. The empirical heat transfer 
correlation of the porosity model needs to be verified in the future. 
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